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(— i 1 Fractional derivative in time variable is introduced into the Fokker-Planck equa- 

^ tion of a population growth model. It's solution, the KNO scaling function, is 

' transformed into the generating function for the multiplicity distribution. Formu- 

las of the factorial moment and the Hj moment are derived from the generating 
function, which reduces to that of the negative binomial distribution (NBD), if the 
fractional derivative is replaced to the ordinary one. In our approach, oscillation of 
Hj moment appears contrary to the case of the NBD. Calculated Hj moments are 
' compared with those given from the data in pp collisions and in e + e~ collisions. 

> : 

ON ■ 

1 Introduction 

o : 

The negative binomial distribution (NBD) is often used for the analysis of 
I , observed multiplicity distributions in high energy hadron-hadron (hh), and 

e + e~ collisions. The Hj moment derived from the generating function of 
the NBD does not show oscillatory behavior as the rank of cumulant moment 
increases. On the othe hand, Hj moments obtained from observed multiplicity 
I | distributions in hh and e + e~ collisions show oscillatory behaviors 0. If the 

truncated NBD is used to calculate the Hj moments in ^.collisions, calculated 
<-j ', results can describe those obtained from the data well&tl However, in e + e _ 

collisions, calculated Hj moments by the use of the truncated NBD oscillate 
. £^ | much weaker than those obtained from the data, and cannot n explain the 

■ oscillatory behavior of the Hj moments obtained from the data □. 

The NBD is derived from a branching process, birth and death process 
with immigration. The probability distribution (the multiplicity distribution) 
and the branching equation are transformed respectively into the probability 
density (the KNO scaling function) and the Fokker-Planck equation by the 
Poisson transform. n 

The fractional calculus has a long history u. Recently, the fractional 
Fokker-Planckcquation have been introduced to describe anomalous diffusion 
phenomena &Lra; the fractional space or time derivative replaces the ordinary 
space or time variable. If fractional derivative in time variable is introduced to 
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the Fokker-Planck equation, it means that a sort of time lag or a memory effect 
is introduced into the Fokker-Planck equation, or the branching equation. 
We would like to investigate this type of fractional Fokker-Planck equation, 
the solution of which reduces to the gamma distribution when the fractional 
derivative is replaced to the ordinary one, and an effect of fractional derivative 
on the behavior of cumulant moments. 



2 The Fractional Fokker-Planck Equation 

In the following we investigate the fractional Fokker-Planck equation (FFPE), 



dt 

Lpp 



= T>\- a L FP ij){z, i), < a < 1 , 
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a(z) = (3 — jz, b(z) = cr 2 z, (1) 



with the initial condition, 

4>(z,t = 0) = 8(z-Zo). (2) 

In Eq.(|), P > 0, a > and jy is real. Symbol oD* denotes the Riemann- 
Liouville fractional derivative El defined by 



oD?/(t) 



1 



d" 
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/(r)dr, n - 1 < S < n, (3) 



T(n - 5) dt n J 

where n is a natural number. If a = 1, oD®f(t) = f(t) ( oD? = 1 ), and 
Eq.(|l]) reduces to the ordinary Fokker-Planck equation. ~ 

According to the method proposed by Barkai and Silbey 19, we assume 

that 



il>{z,t) = / R s {t)G s (z)ds, (4) 
Jo 

and that function G s {z) satisfies the following equations, 

L FP G s (z) = ^G a (z), G (z)=S(z-z ). (5) 

Then function G s (z) becomes the KNO scaling function for the Pefina-McGill 
formula, 



Gs{z) = ^-exp ] 

kp kp 



k 
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Applying the Laplace transform to Eq. ([!]), we find 



uR s (u) + u 



ds 



G s {z)ds = 1 - v}- a Ro{u) S(z - zq), (7) 



where R s (u) is the Laplace transform of R s (t), 

/>00 

R s {u) = / R s {t)e- Ut dt. 
Jo 

Furthermore, we assume that each side of Eq.(^) is equal to zero; 

1 _ a dR s (u) 



v}- a R (u) = 1, 



du 



uR s (u). 



(8) 



(9) 



The solution of Eq.(^) is given by R s (u) — u a 1 exp[— su a }. Then R s (t) is 
written as 



Rs(t) = 77- / ^ cxp[a - -a a ]da, (c > 0). (10) 

In the following, we assume that 7 > 0. Equation (j^j can be expanded as 
G s {z) = I(|)^exp[-f] ^ 7^-TT^- 1) (r)4t- 1) (^)e-" 1 -, (11) 



r(m + A) 



where Lm ^ (z) denotes the Laguerre polynomial. Then the solution for the 
FFPE (0) with z = +0 is given by 

A-l 00 ' 

z z, x m 



t .^ r(m+A) 4;-"(j)^(-^), (12, 

TCI— X 7 



where E a (—t) for £ > denotes the Mittag-Leffler function 0, 



E a (-t) 



sin(a7r) 



exp[— (xt) 



l/ai 



It is written in the infinite series as 



x 2 + 2xcos(«7r) + 1 



E a {z) = J2 



jT(on + l) 



da;. (13) 



(14) 



if a = 1, Eq.(|l^) reduces to the gamma distribution, the KNO scaling function 
of the NBD. 
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3 Factorial Moment and Cumulant Moment 



The generating function for the multiplicity distribution P(n,t) is defined as 

oo 

n(u) = ^p(n,ty\ (15) 

The multiplicity distribution P(n, t) is connected to it's KNO scaling function 
ip(z,t) by the Poisson transform, 

p(n }t ) = ^}l f z n ex.p[-(n )zU(z,t)dz. (16) 
nl Jo 

Then the generating function (GF) is written as 

/•OO 

n(l - u/(no» = / ip(z,t)e- uz dz. (17) 
Jo 

The GF corresponding to Eq.fO) is given by 



*(«) = E r( ?w^ ( ' fc( " 0)(M ' 1))? " g a (-m 7 t"). (18) 
^ m!r(A) (l-fc(no)(n-l)r +A 

The factorial moment are defined as 
dm{u) 



Si 



= <n(n-l)---(n-i + l)). (19) 

U — 1 



<9uJ 

The normalized factorial moment is given from Eqs.© and @ as 

p _ _ r(A + j) S 3 m= o(-l) m jCm^a(-m 7 t°) 

J (n)i r(A)AJ [1-£ Q (-7* Q )F ' 1 J 

The fcth rank cumulant moment is defined by the following equation, 

Kj =d j lnll(u)/du j \ u = 1 . (21) 
Then we obtain the recurrence equation for the Hj moment, 

i ~ 1 F 
Hi = 1, Hj = 1 — 2, j-iCm-i J — m H m , (22) 

— i ^7 

m— 1 J 

where i/j = Kj/fj- 

If a = 1, Eq. (|l7j) reduces to the GF for the NBD, with mean multiplicity 

(n) =fcA(n )(l-e-T t ), 

n NB (u) = [l-«n)/A)(u-l)]-\ (23) 
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The normalized factorial moment and the Hj moment for the NBD are given 
from Eq.(|23|) respectively as, 

_ r(A + j) _ r(A + i)(j-i)! 

^ NB ' J -T(A)Ar' HNB ' j ~ t(\+j) • (24) 



As can be seen from Eqs.(pO[) and (|24|), difference between the normalized 
factorial moment derived from the FFPE (0 < a < 1) and that of the NBD 
(q = 1) is given by Mittag-Leffler functions. 

4 Calculated Results and Discussions 

In order to calculate Hj moments, we use observed values of (n) and C 2 (= 
(?i 2 ) / {n) 2 ) for the charged particles. Then, if a and 7i Q are given, A in Eq.(po|) 
is determined by the following equation, 

1 _ ( C 2 _ ±\ [l-£q(- 7 * a )] 2 _ L 



A V (n) J 1 - 2E a (-jt a ) + E a (-2-/t a ) 

In the following, calculated results with a = 0.5 are shown. In Fig. la, our 
results with ^t a = 1.0 (A =66.82) and 2.0 (10.24) are compared with the H 3 
moment obtained from the data in pp collisions at yfs = 546GeV □. We can 
see that calculated values with a = 0.5 and jt a = 1.0 oscillate as strong as 
those from the data. If parameter ^yt" is adjusted with a step of 0.01 so that 
the first relative minimum of the calculated Hj moment should be located at 
j > 5, the calculated first relative minimum value is Hr — —3.323 x 10 -5 , and 
the absolute value of it is much smaller than that obtained from the data. 

In Fig. lb, the calculated Hj moment is compared, with that obtained 
from the data in e + e~ collisions at y/s — 91 GeV I— I. Parameter jt a is 
adjusted so that the first relative minimum should be located at j > 5. As 
can be seen from the figure, the strength of the oscillation of calculated Hj 
moment is comparable with the data. Even if the truncated NBD is used 
to calculate the Hj moments,_the oscillation of calculated moments is much 
weaker than that of the data □. Therefore, we can improve the results much 
better by introducing fractional derivative into time variable of the Fokker- 
Planck equation. 

Calculated results from the GF with a =0.5 show that the oscillation 
of Hj moments are comparable with the data in amplitude both in pp and 
e + e~ collisions. However, the first relative minimum points, or periods of 
oscillations are not necessarily the same as the data. j—. 

The fractional moments have been also investigated by HegyixjJ However, 
the formula derived from Eq.(hJ) is different from those in ref. Ej. 
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Figure 1. Calculated Hj moments as a function of rank j are compared with Hj moments 
obtained from the data, a) in pp collisions, and b) in e + e — collisions. 
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